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ON  THE  EFFECT  CF  THICKNESS  ON  THE  DAMPING  IN  ROLL 
OF  AIRFOILS  AT  SUPERSONIC  SPEEDS 


ABSTRACT 

The  effect  of  thidmess  on  the  damping  in  roll  of  airfoils  at 
supersonic  speeds  is  i castigated.  The  C and  the  lifting  pressure 

1 P 

distribution  for  several  simple  airfoils  are  found  by  the  use  of  a second 
order  theory  similar  to  the  second  order  theory  introduced  by  Busemann 
and  extended  by  Van  Dyk» . 

The  airfoils  considered  consist  of  infinite  wings  with  arbitrary 
symmetrical  cross  sections  both  s'-’ept  and  unswept,  and  of  a reversed 
delta  wing  which  las  a paraboUc  cross  section. 

The  flows  considered  are  three-dimensional  flows,  and  the  methods 
employed  are  somewhat  different  from  those  used  by  Van  Dyke.  The  possi- 
bility of  using  the  methods  e^loyed  in  the  present  paper  for  other  flow 
problems  is  discussed  briefly. 
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SYMBOLS 


A aspect  ratio 

a local  velocity  of  sound 

b airfoil  span 

c velocity  of  sound  in  free  stream 

c velocity  of  sound  in  air  which  has  been  brougnt  to  rest  adia- 

° batically 

C_  pressure  coefficient  (pE225ure.  ) 

p / _2 

1/2  PV2 


rolling  moment  \ 

1/2  pT2  ^7-sb/ 

d airfoil  chord 

f arbitrary  function  associated  with  the  equation  of  the  airfoil 

surface'  (see  equation  (59)) 

k,  k^,  kg,  k^  constants 

K(t)  arbitrary  function  of  time 

l a distance  small  compared  to  unity 

M Mach  number 

m slope  of  leading  edge  (cot  .A  ) 

N - (y  + 1)  M2 

2P2 

\ - -IP2  \ ♦ >2  ♦ kv3 

P pressure 

p rate  of  roll 

q velocity  vector 


\ 


damping  in  roll  lim 
P o 


lm  / 

-fc-  o ' 


U 


/!/ 


5 )2  - p2  (7 )2  - P2 


r,  e 
s 
8 
t 
V 

X»  7,  * 


p 


polar  coordinates 

area  of  airfoil 

equation  of  airfoil  surface 

tine 

free  stream  velocity 
rectangular  coordinates 

/ p2»2-i 

/7T 


P if  P 2,  P3,  I™*  U auxil lary  functions  used  in  finding  the  second 

order  potential  functions 


T 

£ . 

£ *?•  ? 
0,  0„,  and  6. 


’1  w2» 


.A 

X 


V v2*  v3 


ff7- 

$ 

* 


adiabatic  exponent 

parameter  small  compared  to  unity 

rectangular  coordinates 

auxiliary  functions  used  in  finding  the  second  order 
potential  function 

leading  edge  sweep  back  angle 

minus  the  slope  of  the  center  line  of  a reversed  delta 
wing 

direction  cosines  of  outward  normal  to  element  of  area  da 

density 

c af 

£ a* 

. af 

t ay 

second  order  perturbation  potential  function 

first  order  perturbation  potential 

second  order  correction  potential 
scalar  potential  function 


5 


INTRODUCTION 

The  damping  in  roll  coefficient,  C.  , provided  by  fins  and  wings 

is  important  in  the  design  of  certain  types  of  finned  missiles  and  in 
stability  studies  of  aircraft.  Fortunately  for  many  cases  the  effects 
of  missile  bodies  and  aircraft  fuselages  on  the  C can  be  neglected. 

TP 

In  these  cases  the  C can  be  found  from  a study  of  the  isolated  air- 

foils.  The  use  of  the  linearized  theory  of  supersonic  flow  has  permitted 
an  evaluation  of  the  C for  a number  of  different  types  of  wings  and 

lP 

fins  (for  example  see  references  1 through  6).  Expressions  for  the  C. 

P 

found  by  the  linearized  theory  depend  upon  the  Mach  number  and  the  plan- 
form  geometry  and  are  independent  of  the  thickness  distribution.  Tests 
conducted  in  the  free  flight  Aerodynamic  Range  of  the  Ballistic  Research 
Laboratories  at  Aberdeen  Proving  Ground  (reference  8)  indicate,  however, 
that  for  certain  rectangular  fins  the  thickness  distribution  has  an 
appreciable  effect  on  the  C • Thus,  it  appears  that  a study  of  the 

1 P 

second  order  flow  around  rolling  airfoils  might  yield  some  useful  infor- 
mation. 

The  present  paper  presents  the  C for  several  simple  airfoils  cal- 

*P 

culabed  on  the  basis  of  a second  order  theory  similar  to  the  one  intro- 
duced by  Busemann  (reference  9)  and  extended  by  Van  Dyke  (references 
10  and  11).  The  airfoils  considered  herein  consist  of  infinite  wings 
with  arbitrary  symmetrical  cross  sections  both  swept  and  unswept,  arvi  of 
a reversed  delta  wing  which  has  a parabolic  cross  section.  It  is  the 
authoitf*  opinion  that  fran  the  study  of  these  relatively  simple  airfoils 
the  effect  of  thickness  on  raary  ether  airfoils  can  be  estimated. 

The  analysis  contained  in  the  present  paper  utilizes  the  same  type 
of  iteration  as  that  used  by  Van  Dyke  in  reference  (10  and  11).  There 
are,  however,  two  distinct  differences  between  the  analyses  of  Van 
Dyke’s  and  the  present  paper.  Firstly,  the  partial  differential  equa- 
tions differ;  secondly,  the  methods  of  determining  the  second  order 
solution  are  not  the  same. 


ANALYSIS 


Introduction; 

Recent  work  by  Milton  D.  Van  Dyke  (references  10  and  11)  indicates 
that  second  order  solutions  of  the  partial  differential  equation  of 
steady  supersonic  flow  can  be  obtained  by  the  use  of  iterative  methods'. 
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The  partial  differential  equation  considered  here  is  not  the  equation  of 
steady  supersonic  flow.  It  is  quite  similar,  however,  and  we  shall  assume 
that  the  second  order  solution  can  be  obtained  by  iterative  methods. 

It  will  also  be  assumed  that  the  characteristics  sure  the  same  for  the 
first  and  second  order  solutions  (this  assumption  was  made  in  references  10 
and  11) • For  steady  plane  flow  the  second  order  solution  (reference  10) 
found  by  using  an  iterative  method  based  on  the  preceding  assumption  yields 
the  correct  second  order  pressure  of  the  Busemann  second  order  theory. 
Unfortunately,  no  such  Justification  of 'the  preceding  assumption  is  known  to 
the  authors  for  three  dimensional  flows*. 

The  iterative  method  used  in  references  10  and  11  requires  a parti- 
cular integral  of  the  nonhoroogeneous  second  order  partial  differential 
equation  in  terms  of  the  first  order  solution.  The  particular  integrals 
utilized  in  the  present  paper  apply  only  to  the  particular  problem  under 
consideration.  The  method  developed  here  has  the  advantage  that  it  can 
be  applied  to  problems  in  three  dimensional  flow  where  a particular  inte- 
gral of  the  nonhoraogeneous  second  order  equation  in  terms  of  the  first 
order  solution  is  not  available  in  the  literature. 


The  Partial  Differential  Equation? 

The  partial  different.' rl  equation  for  the  potential  function  to  be 
used  i*i  a special  case  of  the  three  dimensional  time  dependent  equation 
for  the  potential  function  of  a non-viscous  compressible  fluid.  For 
completeness  this  equation  will  be  determined  from  the  equations  of 
hydrodynamics  • 

Four  of  the  equations  needed  are  contained  in  the  equation  of  con- 
tinuity: 

V.  pq  ♦ ||  - 0 (1) 

and  in  the  three  Euler  equations: 

~ VP  ♦ \ V q2  ♦ - o (2) 

p c 9t 

where  q is  the  velocity  vector  and  p is  the  density. 

Equations  (1)  and  (2)  give  four  relations  between  the  five  variables 
P,p  and  the  three  components  of  q.  A fifth  relation  is  obtained  from  the 
assumption  of  *ero  heat  transfer  (constant  entropy).  In  the  case  of  a 
perfect  gas  this  relation  is 

P - k pY  (3) 
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where  k lb  a constant  and  y Is  the  adiabatic  exponent* 

It  will  be  assumed  that  the  velocity  can  be  expressed  as  the  gradient 
of  a potential  function) 


q - V-ft.  (U) 

Equation  (2)  may  therefore  be  written  as 

-vp  + i v q2  ♦ 0 

The  preceding  equation  can  be  integrated  to  yield 

-ftt  * 7 **  * It  ' 1 (t) 

where  K(t)  is  an  arbitrary  function  of  time.  The  function  K(t)  is  a 
constant  (to  be  denoted  by  K)  for  problems  which  will  be  considered  here 
since  the  flow  upstream  of  the  airfoils  Is  independent  of  time* 

The  result  of  evaluating  the  integral  in  the  preceding  equation  by 
the  use  of  equation  (3)  is 

-ft  t *i  pT-l.K  (5) 

4 2 (r-D 

The  first  partial  derivative  of  equation  (5)  with  respect  to  time  can  be 
expressed  as  (with  the  aid  of  equation  (3)  ) 

-n-u*i  §£  ♦?  I - ° <6) 

Equation  (1)  can  be  expressed  as 

• <1+  ^ 1 *VP  ♦ ® (7) 

dP  2 

Denoting  the  square  of  the  velocity  of  sound,  by  a it  follows  that 

equation  (7)  can  be  written  as 

f£  ) - 0 (8) 

Eliminating  P ana  p from  equation  (8)  by  the  use  of  equations  (2)  and 
(6)  yields 

2 

^ * q ■ ♦ q • V q2  (9) 
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Equation  (9)  can  be  written  in  terms  of  the  derivatives  of  the  potential 
function  and  as 

a2  ▼ 2 SL  - iltt  ♦ \ ^ SI  • V ( VJl  • VSl)'  ♦ (vjtt.Vil) 

(10) 

Equation  (10)  is  the  three  dimensional  time  dependent  partial  differential 

2 

equation  for  the  potential  function  provided  a is  expressed  in  terms  of 
the  derivatives  of  the  potential  function. 

By  definition 


By  equation  (3) 

•Z-»TPrl 

2 

Substituting  this  relation  into  equation  (5)  and  solving  for  a yields 
»s  - e/  - *i2y2  .Jl,2  ♦ 2 Jlt)  (11) 

The  function  a is  the  local  velocity  of  sound  and  the  constant  cq  is 

the  velocity  of  sound  in  the  gas  which  has  been  brought  to  rest  adia- 
batically.  It  is  convenient  to  express  a2  in  terns  of  c,  the  velocity 
of  sound  in  the  undisturbed  stream  ahead  of  the-^irtfeil.  It  follows 
from  equation  (11)  that  c is 

c2  - c02  - (*£)  V2 
2 

The  result  of  eliminating  c from  the  preceding  equation  and  equation  (11) 
is  ° 

a2  - c2  (ilx2  *Ii72  + IL2  ♦2iLt-V2)  (12) 

A perturbation  potential  will  be  used  which  is  defined  by 
D.  - V (x  ♦ £ ) 

Note  that  the  perturbation  potential  lias  been  normalised  through  division 
by  the  free  stream  velocity* 


✓ 
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It  will  be  assumed  that  for  the  airfoils  considered  here  the  triple 
products  in  equation  (l£)  can  be  neglected.  Based  on  this  assumption 
equation  (15)  reduces  to 


•fr  & + g ♦ 

xx  77 


$ “ -7  & - -w  $ " ( $ 4 

zz  c xt  c tt  x 


W £ ) f & ♦ <£  ♦ © )+  2 $ $ +2$  $ ♦ 

t'  ' xx  77  yxz/  x xx  7 ay 


(16) 


2 * $i§rt*7  *.  M 

Equation  (16)  can  be  expressed  in  C7lindrical  coordinates  b7 

“P2  $ ♦ § * $ * T7TT2  $ ^ " 7 # 

x*x*  r»r'  r r*  ; 0*9*  c x»t»  <T 


t«t« 


r 


# 


x»x* 


* & 

“ jtfl 


2 $, 


X' 


£ 


x»x« 


0* 


♦ 2 £ 6 

r*  x*r* 


♦ 


nt? 


The  coordinate  axes  of  equation  (17)  are  fixed  in  space.  It  is  more 
convenient,  however,  to  choose  a new  set  of  axes  fixed  to  and  rolling  with 
the  bod7,  the  x-axis  coinciding  with  the  axis  of  roll.  The  relations  between 
the  components  of  the  two  sets  of  axes  are: 

x-x* 

r»r* 

0-e»  - pt* 
t«t* 


n 


In  the  nev  coordinates  equation  (17)  becomes 


sc  rr  r r r 80 


? #tt*  S #9o'"2['T*1)(fx*T  V ^ *J 


; * 


’ (4  ♦ § + | $ ♦A  * \ *‘4  4 ♦ 

\ xx  rr  r r 06/  x xx 


4-  $ *2$_$ 


r 0 *X9  r r xr 


x xfi  Vr  9 99 


. I£ 

V 


$ f * V $ $ * "Tr  $ $ ♦ 7 $ $ 1 

r r6  x ^xt  Vr^  9 0t  v r rtj 


The  flow  over  bodies  rolling  with  constant  velocities  is  independent  of 
timej  thus,  for  the  cases  to  be  considered  here  the  time  derivatives  are 
zero  in  the  preceding  equation.  Terms  such  as 


4 4 

c 88  r 


£ # * «tc  are  of  the  third  order  and  can 

x ^x8 


be  neglected.  It  follows  from  the  preceding  equation  that  the  partial 
differential  equation  (correct  to  the  second  order)  for  the  flow  over 
bodies  with  constant  rates  of  roll  is  given  (in  rectangular  coordinates) 
by 

-p2  § ♦ & * £ - - 3ft  (y  <|>  * s $ ) ♦ 

xx  yy  zz  c xs  xy 


M2  f (y-1)  $ + ^ + 2 & &>  * (18) 

(.  X \ XX  ‘ JJ  zz/  X “ XX 

2 6 * *2  $ $ l 

y xy  z xa  J 
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Equation  (18)  will  be  used  in  the  following  analysis  of  rolling  airfoils# 
Pressure  Relation t 

The  pressure  coefficient  to  be  used,is  a special  case  of  the  pressure 
coefficient  for  a three  dimensional  tine  dependent  flo*% 

The  pressure  coefficient  will  be  defined  tgr 


(t-  - la ) 
\ % <•„) 


The  square  of  the  local  velocity  of  sound  is 


*2 " 

From  this  and  equation  (3)  it  follows  that 

l 

2 dP  . .rl  UY  „ Y 
a r - k r P 


l - k p r*  m £ 


In  the  undisturbed  stream 


P0  .2 
o — c 

po  T 

The  pressure  coefficient  can  now  be  written  as 


2 / 


1 -ai 


From  equations  (20)  the  pressure  can  be  written  as 


_/  «2  \ rSr 


In  the  undisturbed  stream 
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From  the  two  preceding  relatione 


P 

T 

0 *c 


t5r 


Substituting  this  expression  into  equation  (21)  yields 

2 H)  * 

f\ez  / 


C n - 


4 


Substituting  equation  (lU)  into  the  preceding  equations  gives 

*(*v  4 


(22)' 


Equation  (22)  is  the  pressure  coefficient  for  a three  dimensional  unsteady 
flow. 

Expanding  equation  (22)  about  the  free  stream  conditions  yields 


V*2 


u j-  $2-  #2*P2  « \ + .... 

Vxv/yz  x V V* 


In  rolling  coordinates  equation  (23)  becomes 

it 


(23) 


C ■ -2 
P 


^ $ ♦ “7^ ) “ £ *■  " $ 2 ♦ 02  $ 2 ♦ ( 

v x y * ' y * x 


(2U) 


2M2 


# $ * -*T  ^ * **• 

x ' t V*  t 


For  the  flow  associated  with  steady  rolling  equation  (2U)  becomes  (correct 
to  the  second  order) 

Cp“  ‘ 2 ( # +2V  * “ ^ # ) “ $ " # + (25) 

xxyt/ym  x 

Equation  (25)  will  be  used  in  the  following  analysis  of  rolling  airfoils. 
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Method  of  Iteration: 


It  will  be  assumed  that  equation  (18)  can  be  solved  by  an  iteration 
procedure#  The  first  step  is  to  find  the  first  order  solution.  The  first 
order  partial  differential  equation  is  obtained  by  neglecting  the  second 
order  terms  in  equation  (18).  Thus  the  first  order  differential  equation 
is 

-02  <t  ♦ f -0  (26) 

xx  yy  zz 

The  first  order  solution  is  taken  as  the  first  approximation  to  the 
solution  of  equation  (18).  It  is  assumed  that  the  second  approximation 
can  be  found  by  substituting  the  fir3t  order  solution  into  the  rigvt  side 
of  equation  (18)  and  solvirg  the  resulting  nonhomogeneoua  equation,  which 
from  equations  (18)  and  (26'  is 

-gNlr  ♦ 'Xr  ♦ ( j - % 0 \ ♦ 

xx  'yy  zz  C \ x*  xy  J 

(27) 

2M2  |(H-1)  p2j<x  ♦ <ty 

where 

H.iriLid 

2p 

The  solution  of  equation  (27)  will  be  referred  to  as  the  second  order 
solution. 

The  particular  integrals  utilized  in  solving  equation  (27)  will  only 
apply  to  the  problem  under  consideration.  A different  particular  integral 
is  required  for  each  problem. 


Boundary  Conditions? 

Physical  considerations  require  that  the  flow  be  tangent  to  the 
surface,  and  that  all  velocity  perturbations  vanish  upstream  of  the 
airfoil.  These  boundary  conditions  may  be  expressed  mathematically  as: 

<j(  (x,y, z)  • 0 upstream  of  the  airfoil. 


and 


q . X7  s ■ 0 

where  s(x,y,z)  - 0 is  the  equation  of  the  surface  of  the  body. 
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The  equation  of  the  surface  of  the  body  may  also  be  expressed  as 

a - €t  (x,y) 

where  C is  small  compared  with  unity  and  f(x,y)j  thus, 

♦* 

Since  the  velocity^q,  can  be  written  as 

q - i (V  ♦ vgfx  ♦ V)TX)  ♦ j(V5Ty  ♦ Vfy  ♦ pz)  ♦ k (V^2  ♦ vyz  ♦ py) 

it  follows  that  boundary  condition  on  the  body  surface  is  given  by 

-(V  ♦ V0x  ♦ V^x)  0x  - (VJJfy  ♦ ♦ pz)  ay  ♦ V5fa  ♦ ♦ py  - 0 


3f 

where  <xx  ■ € ^ and  ay  “ ^ 15y  * 

The  coordinate  axes  will  be  chosen  so  that  the  airfoil  under  consider- 
ation lies  approximately  in  the  z - 0 plane.  The  boundary  conditions  for 
the  first  order  solution  are  given  by: 

$ (x,y,z)  * 0 upstream  of  the  airfoil. 


and 


Z"0 


(28) 


Similarly  the  boundary  conditions  for  the  second  order  solutiai  are  given 
hy 

- 0 upstream  of  the  airfoil 

and 

ft  I max<tx  I *°y*y  I I <29> 

z«o  Z“0  z-o  z*o 

Unfortunately  for  the  airfoils  considered  in  this  paper  the  first 
order  velocity  components  are  discontinuous  across  the  Mach  sheet  from 
the  leading  edge.  The  effect  of  these  discontinuities  on  the  second 
order  solution  must  be  evaluated. 


1 6 


a Leading  Edge  Mach  Sheet 
With  Zero  Thickness 


z 


b Leading  Edge  Mach  Sheet 
With  Assumed  Thickness 


<K 


x*JY 


c Discontinuity  In  4,  Across 
Leading  Edge  Mach  Sheet 
Along  The  Line  zsz, 


d Plot  Of  ^Across  Leading 
Edge  Mach  Sheet  With 
Assumed  Thickness  Along 
The  Line  z • z. 


FIG.  1.  An  illustration  of  removing  the  discontinuity  in  <t>*  across  the  leading  edge  Mach  sheet 
by  assumli^  the  Mach  sheet  to  have  thickness. 
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In  order  to  evaluate  the  effect  of  the  discow^ar. a.ti«s  in  the  first 
order  velocity  components  we  shall  assume  that  t~.r  ,=_ding  edge  Mach 
sheet  has  a small  but  finite  thickness,  within  vh4  jh  the  first  order 
velocity  components  will  be  made  continuous . Thus,  the  discontinuities 
• through  the  Mach  sheet  are  replaced  by  continuous  functions.  This  process 
is  illustrated  in  figure  1 for  the  velocity  component  in  the  free  stream 
direction  for  an  airfoil  with  leading  edge  perpendicular  to  the  free  stream. 
The  effect  of  the  discontinuities  on  the  second  order  solution  will  be 
found  by  obtaining  the  second  order  solution  within  the  Mach  sheet  and 
then  letting  the  thickness  approach  zero.* 

The  procedure  outlined  in  the  preceding  paragraph  will  be  illustra- 
ted by  considering  the  case  of  an  infinite  rolling  rectangular  airfoil 
having  a constant  initial  slope, 6 » along  its  leading  edge.  Only  the 
Mach  sheet  above  the  airfoil  will  be  considered  since  the  Mach  sheets 
aboTe  and  below  are  of  the  same  form.  The  discontinuities  through  the 
leading  edge  Mach  sheet  depend  only  on  the  initial  slope  of  the  airfoil; 
therefore  the  discontinuities  in  the  first  order  velocity  components 
through  the  Mach  sheet  above  the  airfoil  will  be  the  same  as  those  through 
the  Mach  sheet  above  a flat  rolling  airfoil  having  a constant  slope  . 

The  first  order  potential  function  associated  with  the  flow  over  the 
upper  surface  of  an  infinite  rolling  rectangular  wing  with  a constant 
slope  € is  given  by 

- e ) l DO) 

where  the  airfoil  lies  approximately  in  the  s-0  plane**. 

The  discontinuities  in  the  first  order  velocity  components  through 
the  leading  edge  Mach  sheet  above  the  airfoil  are  from  equation  (30) 

-6) 

Ufy  • 0 

♦ e ) 

It  will  be  assumed  that  the  Mach  sheet  has  a small  thickness  (see 
Figure  1-b).  The  velocity  components  within  the  sheet  will  be  defined  as 


* An  alternate  approach  would  be  to  attach  a small  surface  to  the  leading 
edge  of  the  airfoil  which  would  cause  the  first  order  velocity  components 
to  be  continuous.  The  effect  of  the  discontinuities  would  then  be  evaluated 
by  a limiting  process  in  which  the  width  of  the  small  attached  surface 
would  approach  zero. 

**  Equation  (30)  can  be  found  by  the  integration  of  a source  distribution; 
(see  references  7 and  12)  however,  for  simple  problems  such  as  this  a 
heuristic  method  generally  yields  the  result  with  a minimum  of  effort. 
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_ /*_ r P*  x 

vl /p  ; 


° 

•».  - - (*?  - <■;  ^ 


(31) 


where  / ia  the  thickness  of  the  Mach  sheet  in  the  x direction.  Note  that 
within  the  Mach  sheet 

-b2  a + a * a - 0 

v rxx  r yy . rzz 

and  that  the  velocity  components  are  continuous  functions  in  the  neigh- 
borhood of  it. 


From  equations  (27)  and  (31)  the  differential  equation  within  the 
Mach  sheet  is  found  to  be 


-a2  \jr  * ♦ tJt 

v T xx  T yy  / az 


♦ 


(32) 


It  is  well  known  (references  7 and  13)  that  the  solution  of  the  non- 
homogeneous  equation 

-e2  - F(x,y,«) 

is  given  by 


Y (x>y»*)  • j|  ///  FC&\?z‘ri  ds  d7  d?  * 


1 
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J J ( j da 
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where 


R - ^(x-£  )2  -p2  (7  ->|)2  -P2  (»-9)2, 

_ 2 

n^  » -10  ♦ j Vg  ♦*  kv^  (*i>¥2#*3  are  the  dlxection 

cosines  of  the  outward  normal  to  the  element  of  area  da,  t is  an  arbitrary 
surface  which  cuts  the  forward  Mach  cone  from  the  point  (x,y,a),  v is  the 
volume  inclosed  in  the  forward  Mach  cone  from  the  point  (x,y,s)  and  the 


surface  s. 


and  the  symbol 


denotes  the  finite  part  of  an  integral  as 


defined  by  Hadamard  in  reference  (12).) 


The  result  of  applying  equation  (33)  to  equation  (32)  is 

Y - - £ jjj ir  nr  jfl  ♦ 

T (3U) 


s 


where  the  volume,  v,  is  the  volume  of  the  Mach  sheet  within  the  forward 
Mach  cone  of  the  point  (x,y,z).  The  surface  s will  be  chosen  so  that  it 
is  made  up  of  the  z-0  plane  and  an  arbitrary  surface  upstream  of  the  Mach 
sheet. 

The  integral  over  the  surface  upstream  of  the  forward  Mach  sheet  is 
* jro  so  that  the  surface  integral  in  equation  ( 3U ) can  be  reduced  to 
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The  area  s^  is  the  area  of  the  zmo  plane  which  is  within  the  Mach  sheet 

and  the  forward  Mach  cone  from  the  point  (x,y,z).  From  reference  7 page 
Z z the  preceding  integral  can  be  expressed  as 

f 

sr  //  <s v*  - r? s’  • \ 

8 


It  follows  from  this  expression  that  equation  ( 3U)  can  be  expressed  as 

-6)  + 


f Jf(  J 


n(EV~  - } 
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. 1 f r 

*J  J 


- (4-  P>;  ) f d£  d^  dj  ♦ 


(ss) 


K - ')« 


d&  d? 


In  order  to  evaluate  the  second  order  solution  on  the  down  stream 
side  of  the  Mach  sheet  the  point  (x,y,z)  will  be  located  at  an  arbitrary 
point  on  the  surface  x * p z + f.  The  thickness  of  the  Much  sheet  was 
chosen  to  be  small.  It  follows  that  for  given  value? of  § and  5 the 
variations  in  the  integrand  of  the  volume  and  surface  integrals  are  made 
up  almost  completely  of  the  variation  of  i * This  allows  the  substitution 

E 

of  y in  the  integrands  except,  in  R.  Equation  (35)  can  now  be  reduced 

to 

* P?  ♦/ 
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Integrating  with  respect  to  77  in  preceding  expression  gives 

z P W 

ycp«*/.r. J d5 j {?<?-*>. 
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Performing  the  remaining  integrations  yields 

/ < P*  * /.7,.>  - - - <?  ) , - 


ftv2  Mf-f)5 


(f  - €) 
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Taking  the  limit  a a f approaches  zero  and  replacing  z by  ^ gives 

r . ft  p*  v 

M2  N <f2  x 


X - pz 


2p* 


(36) 


Equation  (36)  is  the  value  of  the  discontinuity  of  the  second  order 
potential  function  associated  with  the  flow  over  the  upper  surface  of  a 
rolling  airfoil  which  has  a leading  edge  with  a constant  initial  slope  and 
zero  aweepback. 


Similarly  the  discontinuity  of  the  second  order  potential  function 
associated  with  the  flow  over  the  upper  surface  of  a rolling  airfoil  which 
has  a sweptback  leading  edge  with  a constant  initial  slope  is  found  to  be 


where 


mx-y-az-0 


^-3a2y2  + 3ayz-z2j 
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_ Np2m2€  2 

2a 


i. 


(37) 
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The  Infinite  Rectangular  Wing  with  an  Arbitrary  Cross  Section: 

In  this  section  the  second  order  potential  function  associated  with 
an  infinite  rectangular  wing  of  arbitrary  cross  section  will  be  determined. 
The  coordinate  axes  are  chosen  so  that  the  y-axis  coincides  with  the  lead- 
ing edge  (see  Figure  2) • The  flow  over  the  upper  and  lower  surfaces 
are  similar;  therefore,  only  the  flow  over  the  upper  surface  will  be 
considered  in  detail.  The  equation  of  the  upper  surface  will  be  represented 
ty 

z - et  (x)  (38) 

where  € is  small  compared  to  unity  and  f (x)  is  of  the  order  of  unity. 


By  inspection 


P,  - (M2  N ♦ p2)f  (x-pz) 

x vp* 

x-p* 

r2  - -tf2  f f(g>  t"  dfi, 


d£> 


The  function  Pj^  was  found  by  assumiiig  a solution  of  the  form 
ru-.  — (N+k)  P 7JE  ♦ x3  ♦ Kg  f (x+p*)3  ♦ 

2pV 


(x-pz)‘ 


] 


where  K^  and  Kg  are  undetermined  constants.  This  leads  directly  to  the 
result 


2 x3 
-V  x ? + 

3P^ 


-i-ny  f(x*pz)3  ♦ (x-pz)3] 

2Up  L J 
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e.  - 2 T,  - &4~  (M^p2  -l)  f (x-p*)  ♦ 
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x-pa 


-«2f  Hz)  f j 
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. £•*&£(-  („p.)3  ♦ (x-p,)3 
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Since  the  second  order  potential  function  is  expressed  by 


the  second  order  potential  function  is  given  by  (from  equations  (39),  (U2), 
and  (U7)). 

§ - $ * 

f,(x.pi).d^!*.Jf.(I-pl)]2  * (U8) 
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For  airfoils  of  symmetrical  thickness  distribution  at  zero  angle 
of  attack  the  pressure  coefficient  on  the  lower  surface  of  the  airfoil 
is  the  same  as  that  on  the  upper  surface  with  p replaced  by  -p.  The 
pressure  difference  coefficient  can  therefore  be  expressed  as 

«-  Cp  “ Cp  (P>  “ Cp  (-P> 

It  follows  from  the  preceding  relation  and  equation  (U9)  that 

Acp  - -^ri+|  A f (x)~!  (50) 

if* 

Note  that  the  second  order  effect  in  equation  (50)  is  a function  only  of 
the  local  slope  of  the  surracd,  iff'  (x)  . 

It  should  be  noted  that  only  terms  multiplied  by  p and  p £ contribute 
to  the  lifting  pressure.  This  is  true  in  general  since  the  nonhomogene ous 
equation  used  in  finding  the  second  order  solution  is  linear  and.  each  term 

2 

on  the 'right  side  of  the  equation  is  multiplied  by  either  p,  £ >p>pfe> 
or  e In  the  following  analysis  terms  not  multiplied  by  p or  pc  will 
therefore  be  neglected. 


A Region  of  Flow  Over  a Reversed  Delta  Wing? 


The  potential  function  associated  with  the  region  of  flow  which  is 
not  affected  by  the  center  section  of  a reversed  delta  wing  will  be  found 
for  a wing  cross  section  given  by 


z ■ <■ 


_/\d  - 2y  -x, 

T x ) 


The  region  of  flow  considered  and  the  cross  section  are  illtistrated  in 
figure  3* 
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a.  Cross  Hotched  Region  Indicates  The  b.  Cross  Section  View  In  Plane  o-o 
Region  Of  Flow  Not  Affected  By 
The  Center  Section 

FIG.  3.  The  reversed  deiU  wing  with  parabolic  cross  section. 
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The  second  ordor  solution  is 

$ ■ 0 * /*•  0 ♦ ♦ ®2  ♦ ®3 

or  from  equations  (52)  and  (56)  (neglecting  terms  not  multiplied  by  p or  pc  ) 


# . ^ (x-M)  - 1 * .<  iyi-Vi!iy  ♦ 
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From  equations  (25)  and  (57)  the  pressure  difference  coefficient  is  given  by 

*cp " ' tfc  (m2n-2)i2  *b0e 

. e gafagL-  - Eg.yftar!  1 (58) 

P^V 
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The  chord-wise  lifting  pressure  is  plotted  in  Figure  (U)  for  a 
reversed  delta  wing  with  a 5 per  cent  thick  center  section  at  a Mach 
number  of  2 for  two  spanwise  stations.  The  chordwise  lifting  pressure 
variations  presented  in  figure  (U)  are  quite  similar  to  those  for 
infinite  wings  with  the  same  cross  sections. 


The  Inf ini Lo  Sweptback  Rectangular  Wing  With  an  Arbitrary  Cross  Section: 

In  this  section  the  effect  of  sweepback  will  be  investigated  by  con- 
sidering the  infinite  sweptback  rectangular  wing  with  an  arbitrary  cross 
seetion.  The  coordinate  axes  are  chosen  as  indicated  in  figure  (5)* 

The  equation  of  the  upper  surface  will  be  represented  by 


* - € f 


(59) 


The  first  order  solution  is 


<jf  • (rax-y-a*) 
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(60) 
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The  second  order  potential  function  must  satisfy  the  nonhctnogeneous 

2 2 

equation  (neglecting  terms  multiplied  by  p and  £ ) 


The  function  Vr  will  be  divided  into  two  parts  and  8g.  The  . 
function  9^  is  taken  to  be  a solution  of  the  nonhcmogeneous  equation, 
therefore  9^  ia  a solution  of  the  homogeneous  equation. 

The  function  9^  was  chosen  to  be 


3U 


Figure  (6)  presents  the  variation  of  the  chord  vise  lifting  pressure  on  a 
wedge  with  a L5°  sweepback  and  with  a 10  per  cent  thickness  at  a Mach 
number  of  / 5 . 


THE  DAMPING  IN  ROLL 

The  pressure  coefficients  found  in  the  preceding  sections  can  be  used 
to  investigate  the  effect  of  thickness  on  the  damping  in  roll  of  some  air- 
foils. In  the  following  sections  some  exact  and  approximate  expressions 
are  presented  for  the  damping  in  roll  of  several  types  of  plan  forms. 


The  Infinite  Rectangular  Wing: 


The  damping  in  roll  coefficient,  C.  , can  be  written  as 

1 P 


C 7 » list 

Lp  p -P-0 


// 

wing 

area 


y A C dx  dy 
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(66) 


For  an  infinite  rectangular  wing 


the  C 


can  be  written  as 


C 

l 


P 


lira 

b—<® 

p — 0 


y A Cp  dx  dy 


(69) 


Substituting  equation  (50)  into  equation  (69)  and  performing  the  indicated 
mathematical  operations  yields 


(70) 
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FIG. 


Perctfrtoge  Of  Distance  Along  Chord 

. The  chordwise  variation  of  the  lifting  pressure  -VAC  /pd,  on  a wedge  with  a 45®  sweep- 
back  and  a 10%  thickness  at  Mach  numbertf>T  p 
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Equation 


(70) 


indicates  that  the 


effect  of  thickness  on  the 


of  an 


infinite  rectangular  airfoil  with  a symmetrical  thickness  distribution 
varies  directly  as  the  thickness  of  the  trailing  edge.  It  can  also  be 
shown  that  to  the  second  order  the  C is  independent  of  the  angle  of 

? 


attack  of  the  airfoil. 


is 


From  equation  (50)  the  center  of  pressure  on  one  panel  (half-wing) 
found  to  be  . 


x • d - 


KM2  f 


f(x)  dx 


(a) 


1 ♦ M2  N 


(71) 


7 - b/3  (b) 

where  x is  measured  from  the  airfoil's  leading  edge  and  y ?*»  vjasured 
from  the  airfoil's  center  line.  For  the  infinite  wedge  thb  preceding 
equations  reduce  to 


x - d/2 
y - b/3 

Figure  (7)  presents  the  variation  of  x/d  with  Mach  number  for  various 
values  of  the  thickness  parameter  for  an  infinite  wing  with  a parabolic 
cross  section. 


The  Reversed  Delta  Win^t 

An  proximate  expression  for  the  damping  in  roll  of  the  reversed 
delta  wing  previously  considered  can  be  obtained  by  assuming  the  expression 
for  the  pressure  distribution  in  the  region  not  affected  by  the  center 
section  i3  valid  over  the  complete  wing. 

Under  this  assumption  the  (approximate) C . of  the  reversed  delta 

lP 

wing  is  given  by  (from  equations  (58)  and  (69)). 
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FIG.  7.  The  variation  of  the  chordwise  location  of  center  of  pressure,  x/d,  with  Mach  number 
for  various  values  of  the  thickness  parameter  for  an  infinite  wing  with  a parabolic 
cross -section. 
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Figure  (8)  presents  the  variation  of  the  C 


with  Mach  number  for 


0 and  5 per  cent  thick  center  sections  for  various  aspect  ratios.  It  can 
be  seen  from  Figure  (8)  that  the  effect  of  thickness  is  quite  small  in 
this  case* 


From  equation  (58)  the  center  of  pressure  is  located  approximately 
at  the  point 
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Figures  (9)  and  (10)  present  the  variation  of  x/d  and  y/b  respectively 
with  Mach  number  for  various  aspect  ratios  for  reversed  delta  wings  with 
0 and  5 per  cent  thick  center  sections. 


The  Infinite  Swept  Wing: 


The  damping  in  roll  of  an  infinite  swept  wing  with  an  arbitrary 
symmetrical  cross  section  can  be  obtained  by  substituting  equation  (67) 
into  equation  (68)  and  performing  the  indicated  mathematical  operations. 
The  C for  this  type  of  airfoil  is  given  by 


2m 

“15 


€ M2  N02a3f(d) 
a?  d 


(7U) 


Equation  (7U)  indicates  that  the  effect  of  thickness  on  the  C.  of  an 

1 P 

infinite  swept  wing  with  an  arbitrary  symmetrical  thickness  distribution 
varies  directly  as  the  thickness  of  the  trailing  edge . 

Figure  (ll)  presents  the  variation  of  the  C,  with  Mach  number  for 


10  per  cent  infinite  wedges  for  various  sweepback  angles. 
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FIG.  9.  The  variation  of  the  chordwise  location  of  center  of  pressure,  x/d,  with  Mach  muni 
for  various  aspect  ratios  for  reversed  delta  wings  of  parabolic  cross  section  with 
and  5%  thick  center  sections. 


Mach  Number 


FIG.  10.  The  variation  of  the  spanwise  location  of  center  of  pressure,  y/b,  with  Mach  number 
for  various  aspect  ratios  for  reversed  delta  wings  of  parabolic  cross  section  with 
zero  and  5 per  cent  thick  center  sections. 
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Mach  Number 

FIG.  11.  The  variation  of  the  damping  in  roll  coefficient,  Cj  , with  Mach  number  for  10%  thick 
infinite  wedges  for  various  sweepback  angles. 
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The  Swept  Tint 


Figure  6 indicates  that  3 panvise  loading  on  an  infinite  swept  wing 
is  not  syntnetrical . An  approximate  expression  for  the  C of  a swept 

tp 

fin  will  therefore  be  presented*  It  is  assumed  that  pressure  distri- 
bution is  the  same  as  that  for  an  infinite  swept  wing  and  that  the  roll 
axis  is  located  as  indicated  in  Figure  12 • 

The  C*  of  a fin  may  be  expressed  as 
lP 


C - -r~ 
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b (y/m)+d 
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o y/m 


yAC  da 
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Substituting  equation  (67)  into  equation  (70)  and  performing  the  indicated 
mathematical  operations  yields 


m2*3  ( NB2n2 

~j~[  It* 


•9/  v 


Figure  13  presents  the  variations  of  the  C with  aspect  ratio 

P 

for  a 10  per  cent  thick  wedge  for  various  sweepback  angles  at  Mach  number 

✓57  This  figure  indicates  that  for  a fin  with  a wedge  cross  section  the 
difference  between  the  C-  of  the  linearized  theory  and  the  C of  the 

6p  <*p 

second  order  theory  increases  with  sweepback  angle. 

Figure  lii  presents  the  variations  of  the  C with  aspect  ratio  for 

p 

a fin  with  a cross  section  at  y ■ 0 (0  per  cent  thick)  given  by 


10d 


for  various  sweepback  angles  at  Mach  number  J This  figure  indicates 
that  for  a fin  with  a parabolic  cross  section  there  is  little  difference 
between  the  values  of  the  linear  theory  C and  those  of  the  second  order 


theory  C 


CONCLUDING  REMARKS 

Most  of  the  airfoils  considered  in  the  preceding  analysis  are  not 
generally  encountered  in  practice;  however,  many  airfoils  do  have  regions 
with  the  same  pressure  distribution  as  airfoils  considered  here.  It  is, 
therefrve,  possible  to  estimate  fairly  accurately  the  effect  cf  thickness 
on  the  C of  many  airfoils  not  considered  in  the  analysis.  As  an  illus- 
lP 

tration  consider  the  C of  a rectangular  airfoil  with  a wedge  cross 

**P 

section  3uch  as  the  rectangular  fins  considered  in  reference  lU,  Figure 
8 of  reference  lU  presents  the  C of  rectangular  fins  with  8 per  cent 


thickness.  An  estimate  of  the  effect  of  thicicness  can  be  made  by  increasing 
the  C of  the  linearized  theory  by  the  same  percentage  as  the  thickness 

increases  the  C for  an  infinite  wing.  Figure  15  presents  the  values 
T*p 

of  the  C ^ 
lP 

values  presented  in  reference  ll*.  It  can  be  seen  from  this  figure  that  the 
agreement  between  theory  and  experiment  has  been  improved  considerably. 


for  an  infinite  wing, 
estimated  in  this  manner  as  compared  to  the  experimental 


The  results  of  the  present  paper  seem  to  indicate  two  general  trends: 


1.  The  effect  of  thickness  an  the  C.  is  quite  small  for  the 

lP 


airfoils  with  zero  thickness  trailing  edges. 
The  effect  of  thj 
sweepback  angle • 


2 o The  effect  of  thickness  on  the  C increases  slightly  with 

*-P 


The  airfoils  considered  in  this  paper  have  jynmetrical  thickness 
distributions..  But  since  the  flows  over  the  upper  and  lower  surfaces 
of  the  airfoils  treated  herein  are  independent  of  each  other,  aerodynamic 
properties  of  airfoils  with  similar  planforms  but  with  unsymmetrical 
thickness  distributions  can  easily  be  determined  from  the  results  contained 
here. 
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The  limitations  of  the  Busemann  sec  cud  order  theory  for  two  dimensional 
flows  have  been  investigated  (see  reference  15).  Since  the  theory  contained 
in  the  present  paper  is  closely  associated  with  the  Busemann  second  order 
theory  it  seems  likely  that  the  results  presented  herein  have  similar 
limitations • 

The  type  of  analysis  used  in  the  present  paper  is  not  restricted  to 
the  study  of  airfoils  with  a constant  rate  of  roll.  An  examination  of  the 
partxal  differential  equations  and  of  the  boundary  conditions  involved 
indicates  that  the  same  type  of  analysis  can  be  applied  to  air fell a 
at  a constant  angle  of  attack,  airfoils  with  a constant  rate  of  pitch, 
and  airfoils  with  simple  unsteady  motions. 


JOHN  C MARTIN 

qU U 


NATHAN  GERBER 
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